##############################################
#### Cliente: TC923 - Julio Cesar Silva    ###
#### In�cio: 20/08/2019                    ###
#### Autora: Luiz Henrique                 ###
##############################################

#### =============================
#### Preparando o R para an�lise
#### =============================

rm(list=ls())
source("Y:/02 - Produ��o Cient�fica - Trabalhos/Fun��es ABG.R")

#### ========
#### Pacotes
#### ========

if(!require(ggplot2)){ install.packages("ggplot2"); require(ggplot2) }
if(!require(psych)){ install.packages("psych"); require(psych) }
if(!require(psy)){ install.packages("psy"); require(psy) }
if(!require(nFactors)){ install.packages("nFactors"); require(nFactors) }
if(!require(Hmisc)){ install.packages("Hmisc"); require(Hmisc) }
if(!require(fmsb)){ install.packages("fmsb"); require(fmsb) }
if(!require(tm)){ install.packages("tm"); require(tm) }
if(!require(tidyr)){ install.packages("tidyr"); require(tidyr) }
if(!require(dplyr)){ install.packages("dplyr"); require(dplyr) }
if(!require(gridExtra)){ install.packages("gridExtra"); require(gridExtra) }

### ========
### Fun��es
### ========

fa1 <- function(x){
  Mat<- cbind(principal(x,1, rotate="varimax")$loadings
              ,principal(x,1, rotate="varimax")$communality
              ,principal(x,1, rotate="varimax")$weights)
  colnames(Mat)<- c("C.F.","Com.","Peso")
  Mat
}

mean2 <- function(x){
  media <- mean(x, na.rm = TRUE)
  media
}

#### ================
#### Banco de dados
#### ================

setwd("Y:/02 - Produ��o Cient�fica - Trabalhos/TC923 - Julio Cesar Silva/Gr�ficos")
Dados.socio <- read.csv2("Y:/02 - Produ��o Cient�fica - Trabalhos/TC923 - Julio Cesar Silva/Bancos de dados/Dados sociodemogr�ficos.csv")
Dados.dom <- read.csv2("Y:/02 - Produ��o Cient�fica - Trabalhos/TC923 - Julio Cesar Silva/Bancos de dados/Dados dom�nios.csv")
Dados.qvfi3 <- read.csv2("Y:/02 - Produ��o Cient�fica - Trabalhos/TC923 - Julio Cesar Silva/Bancos de dados/Dados QVFI-3.csv")
Dados.qvfi7 <- read.csv2("Y:/02 - Produ��o Cient�fica - Trabalhos/TC923 - Julio Cesar Silva/Bancos de dados/Dados QVFI-7.csv")

###====================
### Arrumando vari�veis
###====================

Recursos.0 <- ifelse(Dados.socio$ESCOLA.RECURSOS == 0, "nenhum",
                     ifelse(Dados.socio$ESCOLA.RECURSOS == 1, "acessibilidade",
                            ifelse(Dados.socio$ESCOLA.RECURSOS == 3, "acessibilidade,livrosembraile",
                                   ifelse(Dados.socio$ESCOLA.RECURSOS == 4, "assistente",
                                          ifelse(Dados.socio$ESCOLA.RECURSOS == 5, "acessibilidade,assistente",
                                                 ifelse(Dados.socio$ESCOLA.RECURSOS == 6, "livrosembraile,assistente",
                                                        ifelse(Dados.socio$ESCOLA.RECURSOS == 7, "acessibilidade,livrosembraile,assistente",
                                                               ifelse(Dados.socio$ESCOLA.RECURSOS == 8, "computador",
                                                                      ifelse(Dados.socio$ESCOLA.RECURSOS == 9, "acessibilidade,computador",
                                                                             ifelse(Dados.socio$ESCOLA.RECURSOS == 12, "assistente,computador",
                                                                                    ifelse(Dados.socio$ESCOLA.RECURSOS == 14, "livrosembraile,assistente,computador","acessibilidade,livrosembraile,assistente,computador")))))))))))
Recursos.1 <- Corpus(VectorSource(Recursos.0))
Recursos.2 <- DocumentTermMatrix(Recursos.1)
Recursos.3 <- as.matrix(Recursos.2)
apply(Recursos.3, 2, sum)

acessibilidade.recursos <- Recursos.3[,"acessibilidade,livrosembraile,assistente,computador"]+ 
  Recursos.3[,"acessibilidade,computador"]+
  Recursos.3[,"acessibilidade"]+
  Recursos.3[,"acessibilidade,assistente"]+
  Recursos.3[,"acessibilidade,livrosembraile"]+
  Recursos.3[,"acessibilidade,livrosembraile,assistente"]; Dados.socio$acessibilidade.recursos <- acessibilidade.recursos

assistente.recursos <- Recursos.3[,"acessibilidade,livrosembraile,assistente,computador"]+
  Recursos.3[,"assistente"]+
  Recursos.3[,"livrosembraile,assistente,computador"]+
  Recursos.3[,"acessibilidade,assistente"]+
  Recursos.3[,"livrosembraile,assistente"]+
  Recursos.3[,"assistente,computador"]+
  Recursos.3[,"acessibilidade,livrosembraile,assistente"]; Dados.socio$assistente.recursos <- assistente.recursos

computador.recursos <- Recursos.3[,"acessibilidade,livrosembraile,assistente,computador"]+
  Recursos.3[,"computador"]+
  Recursos.3[,"livrosembraile,assistente,computador"]+
  Recursos.3[,"acessibilidade,computador"]+
  Recursos.3[,"assistente,computador"]; Dados.socio$computador.recursos <- computador.recursos

livrosembraile.recursos <- Recursos.3[,"acessibilidade,livrosembraile,assistente,computador"]+
  Recursos.3[,"livrosembraile,assistente,computador"]+
  Recursos.3[,"livrosembraile,assistente"]+
  Recursos.3[,"acessibilidade,livrosembraile"]+
  Recursos.3[,"acessibilidade,livrosembraile,assistente"]; Dados.socio$livrosembraile.recursos <- livrosembraile.recursos

nenhum.recursos <- Recursos.3[,"nenhum"]; Dados.socio$nenhum.recursos <- nenhum.recursos

Dados.socio$PERIODO.AT..LAZER.2 <- ifelse(Dados.socio$PERIODO.AT..LAZER == 60, 1,
                                          ifelse(Dados.socio$PERIODO.AT..LAZER == 30, 2,
                                                 ifelse(Dados.socio$PERIODO.AT..LAZER == 15, 3, 4)))

Dados.socio$RENDA.FAMILIAR.2 <- ifelse(Dados.socio$RENDA.FAMILIAR == 6, 1,
                                       ifelse(Dados.socio$RENDA.FAMILIAR == 5, 2,
                                              ifelse(Dados.socio$RENDA.FAMILIAR == 4, 3,
                                                     ifelse(Dados.socio$RENDA.FAMILIAR == 3, 4, NA))))

Dados.qvfi3$Questao3 <- ifelse(is.na(Dados.qvfi3$Q3) == TRUE, mean2(Dados.qvfi3$Q3), Dados.qvfi3$Q3)
Dados.qvfi3$Questao9 <- ifelse(is.na(Dados.qvfi3$Q9) == TRUE, mean2(Dados.qvfi3$Q9), Dados.qvfi3$Q9)
Dados.qvfi3$Questao12 <- ifelse(is.na(Dados.qvfi3$Q12) == TRUE, mean2(Dados.qvfi3$Q12), Dados.qvfi3$Q12)
Dados.qvfi3$Questao13 <- ifelse(is.na(Dados.qvfi3$Q13) == TRUE, mean2(Dados.qvfi3$Q13), Dados.qvfi3$Q13)
Dados.qvfi3$Questao15 <- ifelse(is.na(Dados.qvfi3$Q15) == TRUE, mean2(Dados.qvfi3$Q15), Dados.qvfi3$Q15)
Dados.qvfi3$Questao16 <- ifelse(is.na(Dados.qvfi3$Q16) == TRUE, mean2(Dados.qvfi3$Q16), Dados.qvfi3$Q16)
Dados.qvfi3$Questao17 <- ifelse(is.na(Dados.qvfi3$Q17) == TRUE, mean2(Dados.qvfi3$Q17), Dados.qvfi3$Q17)
Dados.qvfi3$Questao18 <- ifelse(is.na(Dados.qvfi3$Q18) == TRUE, mean2(Dados.qvfi3$Q18), Dados.qvfi3$Q18)
Dados.qvfi3$Questao20 <- ifelse(is.na(Dados.qvfi3$Q20) == TRUE, mean2(Dados.qvfi3$Q20), Dados.qvfi3$Q20)
#Dados.qvfi3$Questao21 <- ifelse(is.na(Dados.qvfi3$Q21) == TRUE, mean2(Dados.qvfi3$Q21), Dados.qvfi3$Q21)# Remover essa quest�o, pois s� tem 1 observa��o!
Dados.qvfi3$Questao22 <- ifelse(is.na(Dados.qvfi3$Q22) == TRUE, mean2(Dados.qvfi3$Q22), Dados.qvfi3$Q22)
Dados.qvfi3$Questao23 <- ifelse(is.na(Dados.qvfi3$Q23) == TRUE, mean2(Dados.qvfi3$Q23), Dados.qvfi3$Q23)
Dados.qvfi3$Questao24 <- ifelse(is.na(Dados.qvfi3$Q24) == TRUE, mean2(Dados.qvfi3$Q24), Dados.qvfi3$Q24)
Dados.qvfi3$Questao25 <- ifelse(is.na(Dados.qvfi3$Q25) == TRUE, mean2(Dados.qvfi3$Q25), Dados.qvfi3$Q25)
Dados.qvfi3$Questao26 <- ifelse(is.na(Dados.qvfi3$Q26) == TRUE, mean2(Dados.qvfi3$Q26), Dados.qvfi3$Q26)
Dados.qvfi3$Questao27 <- ifelse(is.na(Dados.qvfi3$Q27) == TRUE, mean2(Dados.qvfi3$Q27), Dados.qvfi3$Q27)
Dados.qvfi3$Questao28 <- ifelse(is.na(Dados.qvfi3$Q28) == TRUE, mean2(Dados.qvfi3$Q28), Dados.qvfi3$Q28)
Dados.qvfi3$Questao29 <- ifelse(is.na(Dados.qvfi3$Q29) == TRUE, mean2(Dados.qvfi3$Q29), Dados.qvfi3$Q29)
Dados.qvfi3$Questao31 <- ifelse(is.na(Dados.qvfi3$Q31) == TRUE, mean2(Dados.qvfi3$Q31), Dados.qvfi3$Q31)
Dados.qvfi3$Questao32 <- ifelse(is.na(Dados.qvfi3$Q32) == TRUE, mean2(Dados.qvfi3$Q32), Dados.qvfi3$Q32)
Dados.qvfi3$Questao33 <- ifelse(is.na(Dados.qvfi3$Q33) == TRUE, mean2(Dados.qvfi3$Q33), Dados.qvfi3$Q33)
Dados.qvfi3$Questao34 <- ifelse(is.na(Dados.qvfi3$Q34) == TRUE, mean2(Dados.qvfi3$Q34), Dados.qvfi3$Q34)
Dados.qvfi3$Questao35 <- ifelse(is.na(Dados.qvfi3$Q35) == TRUE, mean2(Dados.qvfi3$Q35), Dados.qvfi3$Q35)

Dados.qvfi7$Questao3 <- ifelse(is.na(Dados.qvfi7$Q3) == TRUE, mean2(Dados.qvfi7$Q3), Dados.qvfi7$Q3)
Dados.qvfi7$Questao5 <- ifelse(Dados.qvfi7$Q5 == "75/50", 62.5, as.numeric(as.character(Dados.qvfi7$Q5)))
Dados.qvfi7$Questao16 <- ifelse(is.na(Dados.qvfi7$Q16) == TRUE, mean2(Dados.qvfi7$Q16), Dados.qvfi7$Q16)
Dados.qvfi7$Questao17 <- ifelse(is.na(Dados.qvfi7$Q17) == TRUE, mean2(Dados.qvfi7$Q17), Dados.qvfi7$Q17)
Dados.qvfi7$Questao18 <- ifelse(is.na(Dados.qvfi7$Q18) == TRUE, mean2(Dados.qvfi7$Q18), Dados.qvfi7$Q18)
Dados.qvfi7$Questao19 <- ifelse(is.na(Dados.qvfi7$Q19) == TRUE, mean2(Dados.qvfi7$Q19), Dados.qvfi7$Q19)
Dados.qvfi7$Questao20 <- ifelse(is.na(Dados.qvfi7$Q20) == TRUE, mean2(Dados.qvfi7$Q20), Dados.qvfi7$Q20)
Dados.qvfi7$Questao21 <- ifelse(is.na(Dados.qvfi7$Q21) == TRUE, mean2(Dados.qvfi7$Q21), Dados.qvfi7$Q21)
Dados.qvfi7$Questao22 <- ifelse(is.na(Dados.qvfi7$Q22) == TRUE, mean2(Dados.qvfi7$Q22), Dados.qvfi7$Q22)
Dados.qvfi7$Questao23 <- ifelse(is.na(Dados.qvfi7$Q23) == TRUE, mean2(Dados.qvfi7$Q23), Dados.qvfi7$Q23)
Dados.qvfi7$Questao24 <- ifelse(is.na(Dados.qvfi7$Q24) == TRUE, mean2(Dados.qvfi7$Q24), Dados.qvfi7$Q24)
Dados.qvfi7$Questao25 <- ifelse(is.na(Dados.qvfi7$Q25) == TRUE, mean2(Dados.qvfi7$Q25), Dados.qvfi7$Q25)
Dados.qvfi7$Questao26 <- ifelse(is.na(Dados.qvfi7$Q26) == TRUE, mean2(Dados.qvfi7$Q26), Dados.qvfi7$Q26)
Dados.qvfi7$Questao27 <- ifelse(is.na(Dados.qvfi7$Q27) == TRUE, mean2(Dados.qvfi7$Q27), Dados.qvfi7$Q27)
Dados.qvfi7$Questao29 <- ifelse(is.na(Dados.qvfi7$Q29) == TRUE, mean2(Dados.qvfi7$Q29), Dados.qvfi7$Q29)
Dados.qvfi7$Questao30 <- ifelse(is.na(Dados.qvfi7$Q30) == TRUE, mean2(Dados.qvfi7$Q30), Dados.qvfi7$Q30)
Dados.qvfi7$Questao32 <- ifelse(is.na(Dados.qvfi7$Q32) == TRUE, mean2(Dados.qvfi7$Q32), Dados.qvfi7$Q32)
Dados.qvfi7$Questao33 <- ifelse(is.na(Dados.qvfi7$Q33) == TRUE, mean2(Dados.qvfi7$Q33), Dados.qvfi7$Q33)
Dados.qvfi7$Questao34 <- ifelse(is.na(Dados.qvfi7$Q34) == TRUE, mean2(Dados.qvfi7$Q34), Dados.qvfi7$Q34)
Dados.qvfi7$Questao36 <- ifelse(is.na(Dados.qvfi7$Q36) == TRUE, mean2(Dados.qvfi7$Q36), Dados.qvfi7$Q36)
Dados.qvfi7$Questao37 <- ifelse(is.na(Dados.qvfi7$Q37) == TRUE, mean2(Dados.qvfi7$Q37), Dados.qvfi7$Q37)
Dados.qvfi7$Questao38 <- ifelse(is.na(Dados.qvfi7$Q38) == TRUE, mean2(Dados.qvfi7$Q38), Dados.qvfi7$Q38)
Dados.qvfi7$Questao39 <- ifelse(is.na(Dados.qvfi7$Q39) == TRUE, mean2(Dados.qvfi7$Q39), Dados.qvfi7$Q39)
Dados.qvfi7$Questao40 <- ifelse(is.na(Dados.qvfi7$Q40) == TRUE, mean2(Dados.qvfi7$Q40), Dados.qvfi7$Q40)

###============
### Constructos
###============

AP <- data.frame(Q56 = Dados.dom$Q56, Q69 = Dados.dom$Q69, Q70 = Dados.dom$Q70, Q72 = Dados.dom$Q72, 
                 Q73 = Dados.dom$Q73, Q76 = Dados.dom$Q76, Q80 = Dados.dom$Q80, Q85 = Dados.dom$Q85, 
                 Q91 = Dados.dom$Q91, Q92 = Dados.dom$Q92, Q93 = Dados.dom$Q93, Q100 = Dados.dom$Q100)

IPD <- data.frame(Q10 = Dados.dom$Q10, Q12 = Dados.dom$Q12, Q13 = Dados.dom$Q13, Q14 = Dados.dom$Q14, 
                  Q18 = Dados.dom$Q18, Q23 = Dados.dom$Q23, Q24 = Dados.dom$Q24, Q26 = Dados.dom$Q26, 
                  Q38 = Dados.dom$Q38, Q58 = Dados.dom$Q58, Q63 = Dados.dom$Q63, Q64 = Dados.dom$Q64)

CD <- data.frame(Q17 = Dados.dom$Q17, Q19 = Dados.dom$Q19, Q20 = Dados.dom$Q20, Q25 = Dados.dom$Q25, 
                 Q32 = Dados.dom$Q32, Q34 = Dados.dom$Q34, Q36 = Dados.dom$Q36, Q41 = Dados.dom$Q41, 
                 Q42 = Dados.dom$Q42, Q44 = Dados.dom$Q44, Q47 = Dados.dom$Q47, Q50 = Dados.dom$Q50)

Constructos <- cbind(AP, IPD, CD)

###===================
### An�lise descritiva
###===================

set.seed(30)
t1 <- t(apply(AP, 2, basic.np))

set.seed(30)
t2 <- t(apply(IPD, 2, basic.np))

set.seed(30)
t3 <- t(apply(CD, 2, basic.np))

set.seed(30)
ta <- t(apply(Constructos, 2, basic.np))

a.1 <- rep(1, 12)
a.2 <- rep(2, 12)
a.3 <- rep(3, 12)

Grupos <- c(a.1, a.2, a.3)

ta.1 <- ta[,c(2,4,5)]
tb.1 <- data.frame(Itens=factor(rownames(ta.1), levels = rownames(ta.1)), ta.1, Grupos)

plot1 <- ggplot(tb.1, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens), fill = Grupos)) + 
  geom_bar(stat="identity", color="black", position=position_dodge()) + 
  coord_flip() +
  ylim(0, 5)+
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Instrumento reduzido") +
  theme(plot.title = element_text(hjust = 0.5))+ 
  scale_fill_continuous(name = "Constructos", labels = c("Ang�stia parental", "", 
                                                         "Intera��o projenitor-filho disfuncional", "", 
                                                         "Crian�a dif�cil"))

svg("Solic_Plot1.svg", height = 6, width = 8)
plot1
dev.off()

###=================
### An�lise fatorial
###=================

rbind(
fa1(AP)
,fa1(IPD)
,fa1(CD))

rbind(
Val(AP)
,Val(IPD)
,Val(CD))

###=================================================================
### An�lise descritiva das somas e percentis - Indicadores reduzidos
###=================================================================

Indicadores1 <- cbind(SOMA.AP = Dados.dom$SOMA.AP, SOMA.IPD = Dados.dom$SOMA.IPD,
                     SOMA.CD = Dados.dom$SOMA.CD)

Indicadores2 <- cbind(PERCENTIL.AP = Dados.dom$PERCENTIL.AP, 
                      PERCENTIL.IPD = Dados.dom$PERCENTIL.IPD, 
                      PERCENTIL.CD = Dados.dom$PERCENTIL.CD)

set.seed(30)
desc1 <- t(apply(Indicadores1, 2, basic.np))

set.seed(30)
desc2 <- t(apply(Indicadores2, 2, basic.np))

rbind(desc1, desc2)

###==============
### Gr�fico radar
###==============

### Soma

svg("Solicv4_Radar1.svg", height = 8, width = 8)
Data.1 <- c(apply(Indicadores1, 2, mean))
Data.1 <- data.frame(rbind(rep(50,3) , rep(10,3) , Data.1))
colnames(Data.1) <- c("Ang�stia parental", 
                      "Intera��o progenitor-filho\n disfuncional", 
                      "Crian�a dif�cil")

colors_border <- c(rgb(0.243, 0.439, 0.678, 0.8))

colors_in <- c(rgb(0.243, 0.439, 0.678, 0.2))

radarchart(Data.1, axistype=1 , title = "Soma",
           pcol=colors_border , pfcol=colors_in , plwd=4 , plty=1,
           cglcol="grey", cglty=1, axislabcol="grey", caxislabels=seq(10, 50, 10), cglwd=0.8,
           vlcex=0.95, seg = 4)
dev.off()

### Soma

svg("Solicv4_Radar2.svg", height = 8, width = 8)
Data.2 <- c(apply(Indicadores2, 2, mean))
Data.2 <- data.frame(rbind(rep(100,3) , rep(20,3) , Data.2))
colnames(Data.2) <- c("Ang�stia parental", 
                      "Intera��o progenitor-filho\n disfuncional", 
                      "Crian�a dif�cil")

colors_border <- c(rgb(0.243, 0.439, 0.678, 0.8))

colors_in <- c(rgb(0.243, 0.439, 0.678, 0.2))

radarchart(Data.2, axistype=1 , title = "Percentil",
           pcol=colors_border , pfcol=colors_in , plwd=4 , plty=1,
           cglcol="grey", cglty=1, axislabcol="grey", caxislabels=seq(20, 100, 20), cglwd=0.8,
           vlcex=0.95, seg = 4)
dev.off()


###============================================================
### Compara��o dos indicadores com as vari�veis socioecon�micas
###============================================================

### Ang�stia parental

rbind(
  whitney.abg(APind, Dados.socio$GENERO.FILHO.A.)
  ,whitney.abg(APind, Dados.socio$FILHO.ESTUDA.)
  ,whitney.abg(APind, Dados.socio$ESCOLA.ESPECIAL.)
  ,whitney.abg(APind, Dados.socio$acessibilidade.recursos)
  ,whitney.abg(APind, Dados.socio$assistente.recursos)
  ,whitney.abg(APind, Dados.socio$computador.recursos)
  ,whitney.abg(APind, Dados.socio$livrosembraile.recursos)
  ,whitney.abg(APind, Dados.socio$nenhum.recursos)
  ,whitney.abg(APind, Dados.socio$ATIVIDADE.F�SICA.)
  ,whitney.abg(APind, Dados.socio$ATIVIDADE.DE.LAZER)
  ,whitney.abg(APind, Dados.socio$TRABALHA.FORA.)
  ,kruskal.abg(APind, Dados.socio$ESTADO.CIVIL)$`tabela`
  ,whitney.abg(APind, Dados.socio$VIVE.COM.PAI.)
  ,whitney.abg(APind, Dados.socio$TRABALHO.INFORMAL.)
  ,whitney.abg(APind, Dados.socio$BENEFICIO.SOCIAL.)
  ,whitney.abg(APind, Dados.socio$PLANO.SA�DE)
  ,whitney.abg(APind, Dados.socio$GANHOU.PESO.)
  ,whitney.abg(APind, Dados.socio$DOR.DE.CABE�A)
  ,whitney.abg(APind, Dados.socio$TONTURA)
  ,whitney.abg(APind, Dados.socio$N�USEA)
  ,whitney.abg(APind, Dados.socio$INTOLER�NCIA.�.LUZ)
  ,whitney.abg(APind, Dados.socio$TRISTEZA)
  ,whitney.abg(APind, Dados.socio$ESQUECIMENTO)
  ,whitney.abg(APind, Dados.socio$CANSA�O)
  ,whitney.abg(APind, Dados.socio$DORES.NO.CORPO)
  ,whitney.abg(APind, Dados.socio$INS�NIA)
  ,whitney.abg(APind, Dados.socio$DIFICULDADE.PEGAR.NO.SONO)
  ,whitney.abg(APind, Dados.socio$ACORDA.NA.MADRUGADA)
  ,whitney.abg(APind, Dados.socio$QUEDA.CABELO)
  ,whitney.abg(APind, Dados.socio$HIPERTENS�O)
  ,whitney.abg(APind, Dados.socio$AUMENTO.GLICEMIA))

### Intera��o progenitor-filho disfuncional

rbind(
  whitney.abg(IPDind, Dados.socio$GENERO.FILHO.A.)
  ,whitney.abg(IPDind, Dados.socio$FILHO.ESTUDA.)
  ,whitney.abg(IPDind, Dados.socio$ESCOLA.ESPECIAL.)
  ,whitney.abg(IPDind, Dados.socio$acessibilidade.recursos)
  ,whitney.abg(IPDind, Dados.socio$assistente.recursos)
  ,whitney.abg(IPDind, Dados.socio$computador.recursos)
  ,whitney.abg(IPDind, Dados.socio$livrosembraile.recursos)
  ,whitney.abg(IPDind, Dados.socio$nenhum.recursos)
  ,whitney.abg(IPDind, Dados.socio$ATIVIDADE.F�SICA.)
  ,whitney.abg(IPDind, Dados.socio$ATIVIDADE.DE.LAZER)
  ,whitney.abg(IPDind, Dados.socio$TRABALHA.FORA.)
  ,kruskal.abg(IPDind, Dados.socio$ESTADO.CIVIL)$`tabela`
  ,whitney.abg(IPDind, Dados.socio$VIVE.COM.PAI.)
  ,whitney.abg(IPDind, Dados.socio$TRABALHO.INFORMAL.)
  ,whitney.abg(IPDind, Dados.socio$BENEFICIO.SOCIAL.)
  ,whitney.abg(IPDind, Dados.socio$PLANO.SA�DE)
  ,whitney.abg(IPDind, Dados.socio$GANHOU.PESO.)
  ,whitney.abg(IPDind, Dados.socio$DOR.DE.CABE�A)
  ,whitney.abg(IPDind, Dados.socio$TONTURA)
  ,whitney.abg(IPDind, Dados.socio$N�USEA)
  ,whitney.abg(IPDind, Dados.socio$INTOLER�NCIA.�.LUZ)
  ,whitney.abg(IPDind, Dados.socio$TRISTEZA)
  ,whitney.abg(IPDind, Dados.socio$ESQUECIMENTO)
  ,whitney.abg(IPDind, Dados.socio$CANSA�O)
  ,whitney.abg(IPDind, Dados.socio$DORES.NO.CORPO)
  ,whitney.abg(IPDind, Dados.socio$INS�NIA)
  ,whitney.abg(IPDind, Dados.socio$DIFICULDADE.PEGAR.NO.SONO)
  ,whitney.abg(IPDind, Dados.socio$ACORDA.NA.MADRUGADA)
  ,whitney.abg(IPDind, Dados.socio$QUEDA.CABELO)
  ,whitney.abg(IPDind, Dados.socio$HIPERTENS�O)
  ,whitney.abg(IPDind, Dados.socio$AUMENTO.GLICEMIA))

### Crian�a dif�cil

rbind(
  whitney.abg(CDind, Dados.socio$GENERO.FILHO.A.)
  ,whitney.abg(CDind, Dados.socio$FILHO.ESTUDA.)
  ,whitney.abg(CDind, Dados.socio$ESCOLA.ESPECIAL.)
  ,whitney.abg(CDind, Dados.socio$acessibilidade.recursos)
  ,whitney.abg(CDind, Dados.socio$assistente.recursos)
  ,whitney.abg(CDind, Dados.socio$computador.recursos)
  ,whitney.abg(CDind, Dados.socio$livrosembraile.recursos)
  ,whitney.abg(CDind, Dados.socio$nenhum.recursos)
  ,whitney.abg(CDind, Dados.socio$ATIVIDADE.F�SICA.)
  ,whitney.abg(CDind, Dados.socio$ATIVIDADE.DE.LAZER)
  ,whitney.abg(CDind, Dados.socio$TRABALHA.FORA.)
  ,kruskal.abg(CDind, Dados.socio$ESTADO.CIVIL)$`tabela`
  ,whitney.abg(CDind, Dados.socio$VIVE.COM.PAI.)
  ,whitney.abg(CDind, Dados.socio$TRABALHO.INFORMAL.)
  ,whitney.abg(CDind, Dados.socio$BENEFICIO.SOCIAL.)
  ,whitney.abg(CDind, Dados.socio$PLANO.SA�DE)
  ,whitney.abg(CDind, Dados.socio$GANHOU.PESO.)
  ,whitney.abg(CDind, Dados.socio$DOR.DE.CABE�A)
  ,whitney.abg(CDind, Dados.socio$TONTURA)
  ,whitney.abg(CDind, Dados.socio$N�USEA)
  ,whitney.abg(CDind, Dados.socio$INTOLER�NCIA.�.LUZ)
  ,whitney.abg(CDind, Dados.socio$TRISTEZA)
  ,whitney.abg(CDind, Dados.socio$ESQUECIMENTO)
  ,whitney.abg(CDind, Dados.socio$CANSA�O)
  ,whitney.abg(CDind, Dados.socio$DORES.NO.CORPO)
  ,whitney.abg(CDind, Dados.socio$INS�NIA)
  ,whitney.abg(CDind, Dados.socio$DIFICULDADE.PEGAR.NO.SONO)
  ,whitney.abg(CDind, Dados.socio$ACORDA.NA.MADRUGADA)
  ,whitney.abg(CDind, Dados.socio$QUEDA.CABELO)
  ,whitney.abg(CDind, Dados.socio$HIPERTENS�O)
  ,whitney.abg(CDind, Dados.socio$AUMENTO.GLICEMIA))

###======================
### Vari�veis/Indicadores
###======================

### Indicador QVFI3

Saude.geral.da.visao.3 <- cbind(Dados.qvfi3$Q2, Dados.qvfi3$Questao3)
Impacto.familiar.3 <- cbind(Dados.qvfi3$Q4, Dados.qvfi3$Q5, Dados.qvfi3$Q6, Dados.qvfi3$Q7,
                            Dados.qvfi3$Q8, Dados.qvfi3$Questao9, Dados.qvfi3$Q10, Dados.qvfi3$Questao12)
Personalidade.3 <- cbind(Dados.qvfi3$Q11, Dados.qvfi3$Q14, Dados.qvfi3$Questao22, Dados.qvfi3$Questao24,
                         Dados.qvfi3$Questao25)
Compet�ncia.3 <- cbind(Dados.qvfi3$Questao13, Dados.qvfi3$Questao15, Dados.qvfi3$Questao16, Dados.qvfi3$Questao17,
                       Dados.qvfi3$Questao18, Dados.qvfi3$Q19, Dados.qvfi3$Questao20, Dados.qvfi3$Questao23, 
                       Dados.qvfi3$Questao26, Dados.qvfi3$Questao27, Dados.qvfi3$Questao28, Dados.qvfi3$Questao29)
Tratamento.3 <- cbind(Dados.qvfi3$Questao31, Dados.qvfi3$Questao32, Dados.qvfi3$Questao33, Dados.qvfi3$Questao34,
                      Dados.qvfi3$Questao35)


qvfi3.1 <- apply(Saude.geral.da.visao.3, 1, mean2); Dados.qvfi3$qvfi3.1 <- qvfi3.1
qvfi3.2 <- apply(Impacto.familiar.3, 1, mean2); Dados.qvfi3$qvfi3.2 <- qvfi3.2
qvfi3.3 <- apply(Personalidade.3, 1, mean2); Dados.qvfi3$qvfi3.3 <- qvfi3.3
qvfi3.4 <- apply(Compet�ncia.3, 1, mean2); Dados.qvfi3$qvfi3.4 <- qvfi3.4
qvfi3.5 <- apply(Tratamento.3, 1, mean2); Dados.qvfi3$qvfi3.5 <- qvfi3.5

Geral.3 <- cbind(Dados.qvfi3$Q1, qvfi3.1, qvfi3.2, qvfi3.3, qvfi3.4, qvfi3.5)
qvfi3.6 <- apply(Geral.3, 1, mean2); Dados.qvfi3$qvfi3.6 <- qvfi3.6

### Indicador QVFI7

Saude.geral.da.visao.7 <- data.frame(Dados.qvfi7$Q2, Dados.qvfi7$Questao3)
Impacto.familiar.7 <- data.frame(Dados.qvfi7$Q4, Dados.qvfi7$Questao5, Dados.qvfi7$Q6, Dados.qvfi7$Q7,
                                 Dados.qvfi7$Q8, Dados.qvfi7$Q10, Dados.qvfi7$Q11)
Personalidade.7 <- data.frame(Dados.qvfi7$Q9, Dados.qvfi7$Q12, Dados.qvfi7$Q13, Dados.qvfi7$Q14, 
                              Dados.qvfi7$Questao26, Dados.qvfi7$Questao29, Dados.qvfi7$Questao30, Dados.qvfi7$Q31,
                              Dados.qvfi7$Questao32)
Compet�ncia.7 <- data.frame(Dados.qvfi7$Q15, Dados.qvfi7$Questao16, Dados.qvfi7$Questao17, Dados.qvfi7$Questao18,
                            Dados.qvfi7$Questao19, Dados.qvfi7$Questao20, Dados.qvfi7$Questao21, Dados.qvfi7$Questao22,
                            Dados.qvfi7$Questao23, Dados.qvfi7$Questao24, Dados.qvfi7$Questao25, Dados.qvfi7$Questao27, 
                            Dados.qvfi7$Q28, Dados.qvfi7$Questao33, Dados.qvfi7$Questao34)
Tratamento.7 <- data.frame(Dados.qvfi7$Questao36, Dados.qvfi7$Questao37, Dados.qvfi7$Questao38, Dados.qvfi7$Questao39,
                           Dados.qvfi7$Questao40)

qvfi7.1 <- apply(Saude.geral.da.visao.7, 1, mean2); Dados.qvfi7$qvfi7.1 <- qvfi7.1
qvfi7.2 <- apply(Impacto.familiar.7, 1, mean2); Dados.qvfi7$qvfi7.2 <- qvfi7.2
qvfi7.3 <- apply(Personalidade.7, 1, mean2); Dados.qvfi7$qvfi7.3 <- qvfi7.3
qvfi7.4 <- apply(Compet�ncia.7, 1, mean2); Dados.qvfi7$qvfi7.4 <- qvfi7.4
qvfi7.5 <- apply(Tratamento.7, 1, mean2); Dados.qvfi7$qvfi7.5 <- qvfi7.5

Geral.7 <- data.frame(Dados.qvfi7$Q1, qvfi7.1, qvfi7.2, qvfi7.3, qvfi7.4, qvfi7.5)
qvfi7.6 <- apply(Geral.7, 1, mean2); Dados.qvfi7$qvfi7.6 <- qvfi7.6

###================================================================
### Correla��o das vari�veis num�ricas com os indicadores reduzidos
###================================================================

Socio.qvfi3 <- full_join(Dados.socio, Dados.qvfi3, by ="COD")
Socio.qvfi7 <- full_join(Dados.socio, Dados.qvfi7, by ="COD")

DataCor <- cbind(APind, 
                 IPDind,
                 CDind,
                 Dados.socio$IDADE.CRIAN�A,
                 Dados.socio$MINUTOS.AT..FISICA,
                 Dados.socio$FREQ..AT..F�SICA,
                 Dados.socio$PERIODO.AT..LAZER.2,
                 Dados.socio$PESO.DA.M�E..Kg.,
                 Dados.socio$ALTURA.M�E.cm.,
                 Dados.socio$IMC.M�E,
                 Dados.socio$IDADE..M�E,
                 Dados.socio$CARGA.HOR�RIA,
                 Dados.socio$RENDA.FAMILIAR.2,
                 Dados.socio$QTDE.PESO,
                 Dados.socio$TEMPO.PESO,
                 Dados.socio$APETITE,
                 Dados.dom$Soma.LS, 
                 Socio.qvfi3$qvfi3.6, 
                 Socio.qvfi7$qvfi7.6)

rcorr(DataCor, type = "spearman")

###=================================
### An�lise descritiva - 8 a 12 anos
###=================================

Dados.socio812 <- Dados.socio %>% filter(IDADE.CRIAN�A>=8)

### Tabela "2"

rbind(tab(Dados.socio812$GENERO.FILHO.A.),
      tab(Dados.socio812$FILHO.ESTUDA.),
      tab(Dados.socio812$ESCOLA.ESPECIAL.),
      rbind(tab(Dados.socio$acessibilidade.recursos),
            tab(Dados.socio812$assistente.recursos),
            tab(Dados.socio812$computador.recursos),
            tab(Dados.socio812$livrosembraile.recursos),
            tab(Dados.socio812$nenhum.recursos)),
      tab(Dados.socio812$ATIVIDADE.F�SICA.),
      tab(Dados.socio812$ATIVIDADE.DE.LAZER),
      tab(Dados.socio812$PERIODO.AT..LAZER),
      tab(Dados.socio812$TRABALHA.FORA.),
      tab(Dados.socio812$RENDA.FAMILIAR),
      tab(Dados.socio812$ESTADO.CIVIL),
      tab(Dados.socio812$VIVE.COM.PAI.),
      tab(Dados.socio812$TRABALHO.INFORMAL.),
      tab(Dados.socio812$BENEFICIO.SOCIAL.),
      tab(Dados.socio812$PLANO.SA�DE),
      tab(Dados.socio812$GANHOU.PESO.),
      tab(Dados.socio812$APETITE),
      tab(Dados.socio812$DOR.DE.CABE�A),
      tab(Dados.socio812$TONTURA),
      tab(Dados.socio812$N�USEA),
      tab(Dados.socio812$INTOLER�NCIA.�.LUZ),
      tab(Dados.socio812$TRISTEZA),
      tab(Dados.socio812$ESQUECIMENTO),
      tab(Dados.socio812$CANSA�O),
      tab(Dados.socio812$DORES.NO.CORPO),
      tab(Dados.socio812$INS�NIA),
      tab(Dados.socio812$DIFICULDADE.PEGAR.NO.SONO),
      tab(Dados.socio812$ACORDA.NA.MADRUGADA),
      tab(Dados.socio812$QUEDA.CABELO),
      tab(Dados.socio812$HIPERTENS�O),
      tab(Dados.socio812$AUMENTO.GLICEMIA))

### Tabela "3"

rbind(basic(Dados.socio812$IDADE.CRIAN�A),
      basic(Dados.socio812$MINUTOS.AT..FISICA),
      basic(Dados.socio812$FREQ..AT..F�SICA),
      basic(Dados.socio812$PESO.DA.M�E..Kg.),
      basic(Dados.socio812$ALTURA.M�E.cm.),
      basic(Dados.socio812$IMC.M�E),
      basic(Dados.socio812$IDADE..M�E),
      basic(Dados.socio812$CARGA.HOR�RIA),
      basic(Dados.socio812$QTDE.PESO),
      basic(Dados.socio812$TEMPO.PESO))

### Tabela "4"

Dados.dom812 <- Dados.dom %>% filter(IDADE.CRIAN�A>=8)

DI <- cbind(Dados.dom812[,5:13])
RE <- cbind(Dados.dom812[,14:19])
MO <- cbind(Dados.dom812[,20:24])
AC <- cbind(Dados.dom812[,25:31])
AD <- cbind(Dados.dom812[,32:42])
DE <- cbind(Dados.dom812[,43:51])

set.seed(3000)
t1 = t(apply(DI, 2, basic.np))

set.seed(3000)
t2 = t(apply(RE, 2, basic.np))

set.seed(3000)
t3 = t(apply(MO, 2, basic.np))

set.seed(3000)
t4 = t(apply(AC, 2, basic.np))

set.seed(3000)
t5 = t(apply(AD, 2, basic.np))

set.seed(3000)
t6 = t(apply(DE, 2, basic.np))

set.seed(3000)
rbind(t1, basic.np(Dados.dom812$Percentil.DI), 
      t2, basic.np(Dados.dom812$Percentil.RE), 
      t3, basic.np(Dados.dom812$Percentil.MO), 
      t4, basic.np(Dados.dom812$Percentil.AC), 
      t5, basic.np(Dados.dom812$Percentil.AD),
      t6, basic.np(Dados.dom812$Percentil.DE),
      basic.np(Dados.dom812$Percentil.CD))

### Tabela "5"

CO <- cbind(Dados.dom812[,52:64])
AT <- cbind(Dados.dom812[,65:71])
RO <- cbind(Dados.dom812[,72:78])
DP <- cbind(Dados.dom812[,79:87])
SP <- cbind(Dados.dom812[,88:94])
IS <- cbind(Dados.dom812[,95:100])
HE <- cbind(Dados.dom812[,101:105])

set.seed(3000)
t7 = t(apply(CO, 2, basic.np))

set.seed(3000)
t8 = t(apply(AT, 2, basic.np))

set.seed(3000)
t9 = t(apply(RO, 2, basic.np))

set.seed(3000)
t10 = t(apply(DP, 2, basic.np))

set.seed(3000)
t11 = t(apply(SP, 2, basic.np))

set.seed(3000)
t12 = t(apply(IS, 2, basic.np))

set.seed(3000)
t13 = t(apply(HE, 2, basic.np))

set.seed(3000)
rbind(t7, basic.np(Dados.dom812$Percentil.CO),
      t8, basic.np(Dados.dom812$Percentil.AT),
      t9, basic.np(Dados.dom812$Percentil.RO),
      t10, basic.np(Dados.dom812$Percentil.DP),
      t11, basic.np(Dados.dom812$Percentil.SP),
      t12, basic.np(Dados.dom812$Percentil.IS),
      t13, basic.np(Dados.dom812$Percentil.HE),
      basic.np(Dados.dom812$Percentil.PD))

### Tabela "6"

LS <- cbind(Dados.dom812[,106:124])

set.seed(3000)
t14 = t(apply(LS, 2, basic.np))

set.seed(3000)
rbind(t14, basic.np(Dados.dom812$Soma.LS), basic.np(Dados.dom812$Percentil.estresse.vida))

###================
### Gr�ficos barras
###================

t1a = t1[,c(2,4,5)]
t1b = data.frame(Itens=factor(rownames(t1a), levels = rownames(t1a)), t1a)

plot1 = ggplot(t1b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Distrabilidade/Hiperatividade") +
  theme(plot.title = element_text(hjust = 0.5))


t2a = t2[,c(2,4,5)]
t2b = data.frame(Itens=factor(rownames(t2a), levels = rownames(t2a)), t2a)

plot2 = ggplot(t2b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Refor�o parental") +
  theme(plot.title = element_text(hjust = 0.5))


t3a = t3[,c(2,4,5)]
t3b = data.frame(Itens=factor(rownames(t3a), levels = rownames(t3a)), t3a)

plot3 = ggplot(t3b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Humor") +
  theme(plot.title = element_text(hjust = 0.5))


t4a = t4[,c(2,4,5)]
t4b = data.frame(Itens=factor(rownames(t4a), levels = rownames(t4a)), t4a)

plot4 = ggplot(t4b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Aceitabilidade") +
  theme(plot.title = element_text(hjust = 0.5))


t5a = t5[,c(2,4,5)]
t5b = data.frame(Itens=factor(rownames(t5a), levels = rownames(t5a)), t5a)

plot5 = ggplot(t5b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Adaptabilidade") +
  theme(plot.title = element_text(hjust = 0.5))


t6a = t6[,c(2,4,5)]
t6b = data.frame(Itens=factor(rownames(t6a), levels = rownames(t6a)), t6a)

plot6 = ggplot(t6b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) +
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Exig�ncia") +
  theme(plot.title = element_text(hjust = 0.5))


t7a = t7[,c(2,4,5)]
t7b = data.frame(Itens=factor(rownames(t7a), levels = rownames(t7a)), t7a)

plot7 = ggplot(t7b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Compet�ncia") +
  theme(plot.title = element_text(hjust = 0.5))


t8a = t8[,c(2,4,5)]
t8b = data.frame(Itens=factor(rownames(t8a), levels = rownames(t8a)), t8a)

plot8 = ggplot(t8b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Proximidade") +
  theme(plot.title = element_text(hjust = 0.5))


t9a = t9[,c(2,4,5)]
t9b = data.frame(Itens=factor(rownames(t9a), levels = rownames(t9a)), t9a)

plot9 = ggplot(t9b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Restri��o de papel") +
  theme(plot.title = element_text(hjust = 0.5))


t10a = t10[,c(2,4,5)]
t10b = data.frame(Itens=factor(rownames(t10a), levels = rownames(t10a)), t10a)

plot10 = ggplot(t10b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Depress�o") +
  theme(plot.title = element_text(hjust = 0.5))


t11a = t11[,c(2,4,5)]
t11b = data.frame(Itens=factor(rownames(t11a), levels = rownames(t11a)), t11a)

plot11 = ggplot(t11b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Relacionamento com esposo(a)/parceiro(a)") +
  theme(plot.title = element_text(hjust = 0.5))


t12a = t12[,c(2,4,5)]
t12b = data.frame(Itens=factor(rownames(t12a), levels = rownames(t12a)), t12a)

plot12 = ggplot(t12b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Isolamento") +
  theme(plot.title = element_text(hjust = 0.5))


t13a = t13[,c(2,4,5)]
t13b = data.frame(Itens=factor(rownames(t13a), levels = rownames(t13a)), t13a)

plot13 = ggplot(t13b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,5) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Sa�de") +
  theme(plot.title = element_text(hjust = 0.5))

t14a = t14[,c(2,4,5)]
names(t14a) <- c("")
t14b = data.frame(Itens=factor(rownames(t14a), levels = rownames(t14a)), t14a)

plot14 = ggplot(t14b, aes(x=Itens, y=M�dia, scale_x_discrete(limits = Itens))) + 
  geom_bar(stat="identity", color="black", fill="#3462A8", position=position_dodge()) + 
  coord_flip() +
  ylim(0,3) +
  geom_errorbar(aes(ymin=L.I, ymax=L.S), width=.2, position=position_dodge(.9)) +
  labs(title="", x="Itens", y = "Escala") +   
  ggtitle("Estresse da vida") +
  theme(plot.title = element_text(hjust = 0.5))

# Gr�fico 1 - Dom�nio da crian�a completo

svg("Graf1_L.svg", height = 5, width = 10)
grid.arrange(plot1, plot2, ncol=2)
dev.off()

svg("Graf2_L.svg", height = 5, width = 10)
grid.arrange(plot3, plot4, ncol=2)
dev.off()

svg("Graf3_L.svg", height = 5, width = 10)
grid.arrange(plot5, plot6, ncol=2)
dev.off()

# Gr�fico 3 - Dom�nio dos pais completo

svg("Graf4_L.svg", height = 5, width = 10)
grid.arrange(plot7, plot8, ncol=2)
dev.off()

svg("Graf5_L.svg", height = 5, width = 10)
grid.arrange(plot9, plot10, ncol=2)
dev.off()

svg("Graf6_L.svg", height = 5, width = 10)
grid.arrange(plot11, plot12, ncol=2)
dev.off()

svg("Graf7_L.svg", height = 5, width = 5)
plot13
dev.off()

# Gr�fico 5 - Dom�nio do estresse da vida

svg("Graf8_L.svg", height = 5, width = 5)
plot14
dev.off()

###===============
### Gr�ficos radar
###===============

# Gr�fico 2 - Dom�nio da crian�a completo

Dom.crianca.c <- cbind(Dados.dom812$Percentil.DI, 
                       Dados.dom812$Percentil.RE, 
                       Dados.dom812$Percentil.MO, 
                       Dados.dom812$Percentil.AC,
                       Dados.dom812$Percentil.AD,
                       Dados.dom812$Percentil.DE)

svg("Radar1_L.svg", height =7, width=7)
Data.1 <- c(apply(Dom.crianca.c, 2, mean))
Data.1 <- data.frame(rbind(rep(100,6) , rep(0,6) , Data.1))
colnames(Data.1) <- c("Distrabilidade/\nHiperatividade", "Refor�o\nparental", "Humor",
                      "Aceitabilidade", "Adaptabilidade", "Exig�ncia")

colors_border=c(rgb(0.243, 0.439, 0.678, 0.8))

colors_in=c(rgb(0.243, 0.439, 0.678, 0.2))

radarchart(Data.1, axistype=1 , title = "Indicador completo do dom�nio da crian�a",
           pcol=colors_border , pfcol=colors_in , plwd=4 , plty=1,
           cglcol="grey", cglty=1, axislabcol="grey", caxislabels=seq(0,100,25), cglwd=0.8,
           vlcex=0.95, seg = 4)
dev.off()

# Gr�fico 4 - Dom�nio dos pais completo

Dom.pais.c <- cbind(Dados.dom812$Percentil.CO, 
                    Dados.dom812$Percentil.AT, 
                    Dados.dom812$Percentil.RO, 
                    Dados.dom812$Percentil.DP,
                    Dados.dom812$Percentil.SP,
                    Dados.dom812$Percentil.IS,
                    Dados.dom812$Percentil.HE)

svg("Radar2_L.svg", height =7, width=7)
Data.2 <- c(apply(Dom.pais.c, 2, mean))
Data.2 <- data.frame(rbind(rep(100,7) , rep(0,7) , Data.2))
colnames(Data.2) <- c("Compet�ncia", "Proximidade", "Restri��o\nde papel", "Depress�o", 
                      "Relacionamento com\nesposo(a)/parceiro(a)", "Isolamento", "Sa�de")

colors_border=c(rgb(0.243, 0.439, 0.678, 0.8))

colors_in=c(rgb(0.243, 0.439, 0.678, 0.2))

radarchart(Data.2, axistype=1 , title = "Indicador completo do dom�nio dos pais",
           pcol=colors_border , pfcol=colors_in , plwd=4 , plty=1,
           cglcol="grey", cglty=1, axislabcol="grey", caxislabels=seq(0,100,25), cglwd=0.8,
           vlcex=0.95, seg = 4)
dev.off()



















